General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



CR-1 W80 


\ 


(NASA-CB-144S80) A LINEAR SOLOTION OF LIFT N76-27165 

INTERFERENCE IN SQUARE TUNNELS WITH SLOTTED 
TEST SECTIONS OF FINITE LENGTH (General 

Applied Science Labe. , Inc.) 27 p HC $4.00 Onclas 

CSCL 01 A G3/02 4449 



General Applied Science Laboratories inc 



A LINEAR SOLUTION OF LIFT INTERFERENCE 
IN SQUARE TUNNELS WITH SLOTTED TEST 


NASA CR-UA 980 


SECTIONS OF FINITE LENGTH 
By Charles Ruger and Paolo Baronti 


Prepared under Contract No. NASI-13616 by 
GENERAL APPLIED SCIENCE LABORATORIES, INC. 
Westbury, New York 

for 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 



TABLE OF CONTENTS 


INTRODUCTION 
METHOD OF SOLUTION 
NUMERICAL RESULTS 


CONCLUDING REMARKS 



LIST OF ILLUSTRATIONS 


Page 

FIG. 1. SCHEMATIC OF MODEL-TUNNEL CONFIGURATION k 

FIG. 2. LIFT INTERFERENCE PARAMETER ALONG TUNNEL AXIS 

FOR SEVERAL WALL CONFIGURATIONS 8 

FIG. 3. INTERFERENCE PARAMETER ALONG TUNNEL AXIS FOR 

SEVERAL LENGTHS OF SLOTTED SECTION 10 

FIG. k, FLOW VELOCITY DISTRIBUTIONS ALONG UPPER AND 

LOWER WALLS 11 

FIG. 5. FLOW VELOCITY DISTRIBUTION ALONG SIDEWALLS 12 

FIG. 6. EFFECT OF SLOT PARAMETER ON LIFT INTERFERENCE 13 

FIG; 7. EFFECT OF AN INCREASE OF BOTTOM PLENUM PRESSURE 

ON INTERFERENCE PARAMETER 15 

FIG. 8. EFFECT OF INCREASE OF BOTTOM PRESSURE ON OUT- 
FLOW FROM LOWER WALL 16 

FIG. 9. CENTERLINE PRESSURE DISTRIBUTION INDUCED BY 

MODIFICATION OF SOLID WALLS DOWNSTREAM OF SLOTS 1? 

FIG. 10. EFFECT OF MODIFICATION OF SOLID WALLS DOWNSTREAM 

OF SLOTS ON INTERFERENCE PARAMETER l8 



LIST OF SYMBOLS 


Cl 

lift coefficient 

Cp 

pressure coefficient 

c 

tunnel cross sectional area 

• 

h 

tunnel semi -height 

K 

slot parameter 

1 

length of slotted section 

Mop 

free stream Mach number 

P« 

free stream pressure 


plenum pressure 

P 

porosity parameter 

S 

wing span 

u,v,w 

velocity components 

x,y,z 

Cartesian coordinates 

B 

compressibility parameter, (1 - 

r 

vortex strength 

5 

lift interference factor, wC/2r: 

a 

source strength 

$ 

perturbation velocity potential 


length of panels 

Subscript 

1 

wall induced 

m 

model induced 

im 

induced by image vortices 

x,y,z 

denotes partial differentiation 

a' 

induced by wall singularities 

Superscript 

1 

denotes incompressible quant itii 



ABSTRACT 


A numerical program has been developed to determine the effects on lift 
Interference of the finite length of the slotted portion of a square tunnel. 
Also analyzed by the program are the effects produced by small variations of 
static pressure between upper and lower plenum, and by boundary layer develop** 
ment on upper and lower walls downstream of the slotted section^ In the flow 
model, based on compressible linear theory, the wing is represented by a single 
vortex and the wall induced perturbations are obtained by superimposing the 
effects of a series of images of the central vortex and of a distribution of 
sources and sinks on panels along the tunnel walls. The present results can be 
indicative of the magnitude and of the trends of the interference patterns as- 
sociated with slotted test sections of finite length at high subsonic Mach 
numbers. 
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INTRODUCTION 


The classical approach to the study of tunnel wall Interference in 
subsonic flow is to consider test sections of infinite length, represent 
the model by appropriate singularities (vortices, doublets, etc.)* and ob- 
tain closed form solutions for the interference parameters by satisfying the 
wall boundary conditions for the different cases of solid, porous and slotted 
walls. In actual wind tunnels which have ventilated test sections of finite 
length, the interference experienced by the model may be altered to a certain 
degree by the discontinuities in flow condi lions at the beginning and at the 
end of the ventilated portion. This may be particularly true when testing 
lifting three-dimensional models since the large velocity components created 
by the trailing vortices not only produce conditions of inflow and outflow 
between the plenum and the test section but, subsequently, are constrained 
insirie ^he solid portion of the tunnel downstream of the ventilated section. 

In addition, the low momentum air that enters the tunnel from the upper 
plenum and is accelerated by the tunnel flow may decrease the momentum of 
the tunnel flow and create interference pressure gradients. At the same 
cime, the flow of tunnel air into the plenum may lose momentum when mixing 
with the plenum air and thus establish pressure differentials between top 
and bottom portions of theplenum, changing locally, although slightly, the 
plenum pressure with respect to the assumed free stream static value. 

Of current concern is the quantitative determination of these effects 
for the study of three-dimensional flow fields in transonic tunnels. Present 
capabilities in transonic numerical analysis could allow for such a determina- 
tion, e.g.. Reference (1). However, in order to illustrate the important 
trends of the phenomena under consideration and single out the relevance of the 
parameters effecting the interference, simple and more readily computable flow 
models based on a compressible linear theory can be utilized. For this reason, 
linearized solutions have been considered for the flow fields of interest in 
the present investigation. 

Considered herein is the case of a lifting configuration (where an unswept 
wing is represented by a single horseshoe vortex) in a square tunnel with slot- 
ted walls of finite length on all four sides. The effects on lift interference 
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generated by the finite length of the ventilated portion of the test section 
and by small differences between upper and lower plenum pressure are investi- 
gated. Also investigated are the effects produced by variations of wall bound- 
ary conditions downstream of the slots to simulate synwnetric and asymmetric 
(with respect to tunnel centerline) boundary layer development along the upper 
and lower solid walls downstream of the slots. 

The solution technique that has been adopted for the consideration of 
slotted test sections of finite length is a modification of the numerical 
method of Keller and Wright and of Keller^^\ A double series of images 
of the central vortex about the four tunnel walls is first Introduced, as if 
the walls were solid, so as to nearly satisfy the solid wall boundary condi- 
tions upstream and downstream of the slotted region. The tunnel walls are 
then divided into rectangular panels which are represented by source or sink 
distributions of unknown strength. The appropriate boundary conditions are 
satisfied at the centroid of each panel leading to a matrix equation whose 
solution yields the panel source strength. The wall panel strength distribu- 
tions, together with the image vortices, are then used to calculate the wall 
interference at any point in the tunnel. The introduction of image vortices 
permits an improved satisfaction of the wall boundary conditions upstream and 
downstream of the slots with fewer wall panels. 

This numerical method was employed since analytical solutions, e.g. 

References (A) and (5), are limited in general to infinite-length test sections 
and would not readily allow for the combination of ventilated and solid wall 
boundary conditions or for the more general wall boundary conditions considered 
herein. In the present analysis, the surface singularity distribution is taken 
constant in each panel, which is appropriate for the integrated slotted wail 
condition, as established in Reference (2). The model and analysis have the 
flexibility to allow for the inclusion of other wall boundary conditions and 
distributions of wall singularities (as in Reference 3), as well as of a more 
general representation of tunnel models, such as lifting wings and wi ng-centerbody- 
tail conf igurat ions, by surface distributions of elementary vortices and of other 
singularities. 


2 



METHOD OF SOLUTION 


Following standard methods, and In the coordinate system of Figure (1), 
the problem Is expressed In the term of the linearized compressible equation 


XX yy 22 


= 0 


where 


$ $ + $. 
m I 


being the perturbation potential of the flow about the model In free air 
and being the potential induced by the wall constraints. The basic config- 
uration that is exanined Is the one of Figure (1), with the four tunnel vvells 
slotted In the regloi of the model and solid upstream and downstream. The 
boundary conditions in the slotted portion extending fr*'m to assumed 

to be of the form 


i - i <1, < X < 

( 1 ) 

$+K«i> =0. y = + h, )t. <x<Ji, 

where K Is the slot parameter. In the solid region they are 
$ - 0 2 = + h 

( 2 ) 

Modifications to these boundary condi tions, to account for different pressure 
between upper and lower plenum and different shapes of solid walls, are de- 
scribed later. 


The model considered is an unswept f I n I te wing wh I ch Is represented by a 
single horseshoe vortex of strength F Lo picvlde the value of at the tunnel 
walls. The Induced potential •I*. Is provided by a proper distribution of singu 
laritles at the wall or external to it, which satisfy, together with the model 
Induced perturbations, the boundary conditions (l) and (2). 
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FIGURE 1. SCHEMATIC OF MODEL-TUNNEL CONFIGURATION 


The solution is sought by casting the problem into pn Incompressible form 
utilizing the compressibility transformations 

X* * x/23h y' * y/2h z’ «* z/2h 

and with* 

K' « K and « 3F 

where the prime quantities denote the incompressible values. 


The velocity potential induced by the vortex of strength F* and span s' 
is taken as 



city induced by attached and trailing vortices, as given in Reference (6). 


The distribution and strength of the singularities providing the inter- 
ference potential '!>. is obtained as follows. A double series of images of 
the central vortex is first taken as if the walls were solid. Each image vor- 
tex has a perturbation potential of the form of Equation (3) with the 
proper sign and the proper distances. Then, following the procedure of Refer- 
ences (2) and (3), the tunnel walls are divided into longitudinal strips, each 
of which is divided into a number of rectangular panels of arbitrary length. 
Each panel is represented by a source (or sink) distribution of constant but 
unknown strength o'. The potential Induced by each element of strength 
o' is given in Reference (2), The strengths of the vortex images alternate *n 
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sign In the z-direction and are all positive In the y-dlrettlon. Slncj the 
tunnel under consideration consists of solid walls upstream and downstream of 
the slotted section, use of even a finite number of vortex reflections will 
cause the boundary conditions on the panels in the solid region to be nearly 
satisiied. This results in small source strength distributions oi, these 
panels and allows for use of larger and, therefore, fewer panels in the solid 
region. Use of vortex reflections also results in the solid wall boundary con 
dition being nearly satisfied up and downstream of the axial range in which 
panels are taken- 

As discussed in Reference (2), for a constant strength distribution in 
each element, the integrated form of the slotted wall boundary condition is 
appropriate. Accord i- ply integration of Equation (l), yields 


o — z 




23h 


< x' < 




23h 


I „ 


+ K' , ) * 0 

o ~ y y fC 


where subscript o refers to quantities at the beginning of the slotted region. 

Note that the similarity variable describing the finite slotted lengths are 

A^/23h and Thus, satisfaction of Equations (2) and (^) In terms of 

the model potential (known) and the potentials $I (known) and given 

m im a 

by the image vortices and by the wail singularities, respectively, yields a 
set of linear homogeneous equations, which by matrix inversion give the un- 
known panel strengths a’. From these and from the Image vortices, and by ap- 
plication of the compressibility transfer' it ions in an inverse manner, the 

2 

wall induced interference in terms of u = u'/^ , v = v'/B and w = w'/3 is 
computed at any point in the tunnel, l.ie lift interference is customarily 
cast in terms of the parameter 6 = vt.C/2Ts where w. is the induced velocity 
along ‘•he tunnel axis. This parameter Is invariant under transformation. 

Thus, th? axial distribution of 5 in terms of x' = x/23h and of the scaled 
slotted test section lengths i^/23h and (see Equation ^i) is similar 

(and thus Independent specifically of Mdch number and tunnel height). 
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In carrying out the nuirerlcal csxnputatloni , eight reflections of the 
cehtraTVof tex^^ a been used. Also, the panels are chosen 

such that the stations at which the wall changes from closed to slotted are 
located at tne panel boundary. In order to expedite the numerical solution, 
further advantage wi. taken of the symmetries of the problem. For the prob- 
lem cf a horseshoe vortex in a slotted-sol id tunnel with straight walls and 
equal slots on all four sides not only is there a y-synmetry but, in addition, 
the problem is anti-symmetric with respect to z. In fact, the source strength 
distributions on the bottan wall ;jnd on the lower half of the sidewall are 
equal to and of opposite sign from the distributions on the upper well and the 
upper ha’^ of the sidewall, respectively. This is evident when is substi- 
tuted In the boundary conditions (1) and (k) . Under the above conditions it is 
only necessary to compute the wall source strength distributions at penels a- 
long one quarter of the tunnel cross section (i.e. one-half along the top and 
one-half along the sidewall). Thus, the number of equations to solve Is re- 
duced by a factor of four, although the influence of all the panels around 
the tunnel must be in .luded in the boundary conditions at the centroids of the 
panels chosen. For the case where the lower wall plenum pressure is perturbed 
or the downstream upper solid wall is curved, only the y-syiwnetry remains; still, 
tlw wall sources can be computeu in the half plane, that is along one-half of 
the tunnel cross section. 


NUMERICAL RESULTS 

The typical configuration considered here is a lifting wing of unit span 
and lift coefficient *» 0.4 in a square tunnel of height three times the 
span, at a Mach number of 0.5* Shown in Figure (2) is the lift Interference 
parameter, 5, along the axis of the tunnel for completely solid, completely 
slotted and partially slotted wall; . In all cases, eight reflections around 
each wail of the horseshoe vortex representing the wing were used with strips 
extending over the range of x from -3 to 4 (all dimensions are scaled to the 
wing of unit span). The partially slotted tunnel considers slots extending 
from X “ -2 to X « 2 with a slot parameter K » 0.5. These computations have 
been carried out in a quarter-plane. Several combinations of strips and panel 
lengths were used tc indicate the sensitivity of the solution to the number c.^ 
strips and panels. As indicated in the figure a rather small number of strips 
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Solid Walls, 8 Strips, Ax 
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FIGURE 2. LIFT INTERFERENCE PARAMETER ALONG TUNNEL AXIS FOR SEVERAL 
WALL CONFIGURATIONS 





and panels is sufficient to accurately predict wall interference along the 
tunnel axis. Also indicated in the figure are comparisons at x * 0 with 
other known solutions (References 6 and 7) • 

In the case of partially slotted walls, the drastic change in the inter- 
ference pattern in the region near the downstream transition from slotted to 
solid wall was expected. What is of practical importance is the extent of 
this region in terms of wing span and of tunnel height. The effect of short- 
ening the downstream extent of the slotted walls is shown in Figure (3) where 
results are presented in terms of the similarity variables x/23h and l/2^h. 

Not only can the effect of the transition from slotted to solid wall be felt 
at the wing location, if the wing is relatively close to the end of the slot- 
ted section, but more importantly, in the case of a model with a tail, the 
latter may experience a different wall interference than the wing. Thus, the 
pitching moment of the model may be directly affected by the finite length of 
the slotted test section. 

Complimentary to the results of Figure (2) are those presented in Figures 
(^)., (5) and (6). Shown in Figure {k) are axial distributions of the flow 
velocity (in terms of the parameter wC/2rs) through the slotted portion of 
the walls at two selected locations along tiie top and bottom walls. Similar 
distributions on the sidewalls are shown In Figure (5)* Figure (6) shows the 
variation of lift interference for several values of the slot parameter K. 

The result reiterates, as is well known, the sensitivity of 6 to values of K, 

(as well as to the parameters P and n when a more general boundary condition of 
the type ^ = 0 is used, c.f. References 3 and 8). 

As indicated in some detail in Reference (8), mass flow into the plenum 
may alter plenum pressure conditions with respect to the usually assumed free 
stream static pressure value. Within the present formulation it is possible 
to examine the effects of a small perturbation in the plenum pressure. A com- 
putation has been performed by increasing the pressure in the bottom plenum to 

a constant value of 1.002 p (pressures of this order were observed at certain 

00 

plenum locations during the experiment reported in Reference 9), while retain- 
ing a value of p = p^ along the other walls. Under these conditions, the bound- 
ary condition of Equation (l) can be rewritten, along the bottom wall, in the form 
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FIGURE 3. INTERFERENCE PARAMETER ALONG TUNNEL AXIS FOR SEVERAL LENGTHS OF 
SLOTTED SECTION 
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FIGURE 4. FLOW VELOCITY DISTRIBUTIONS ALONG UPPER AND LOWER WALLS 






Slotted Length 
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FIGURE 6, EFFECT OF SLOT PARAMETER ON LIFT INTERFERENCE 



and the computational procedure carried out in the usual manner. The effects 
on lift interference of this increase in pressure are significant, as shown in 
Figure (7), where the comparison between this situation and the one where the 
plenum pressure is taken constant on all four sides is shown. At the same time, 
the outflow through the bottom wall is greatly reduced, as typically shown in 
Figure (8), whereas the upper wall inflow is only slightly reduced. 

Inflow of stagnant air from the upper plenum into the tunnel will also 
cause an asymmetric development (with respect to tunnel centerline) of the vis- 
cous flow layer along the solid walls downstream of the slotted section. To in- 
vestigate, at least in a qualitative manner, this effect and also to investigate 
the magnitude of a centerline axial pressure gradient which may be caused by a 
variation in cross sectional area downstream of the slots, two computations 
have been performed with modified solid walls and with the same lifting wing of 
unit span and C|_ = 0.4 in the tunnel. In one case, the upper downstream solid 
wall was assumed to converge from x = 1.75 (end of the slots) to x = 4 by an 
amount of 0.05, a number consistent with expected boundary layer growth. In 
the other case, both upper and lower walls were taken as convergent. The side- 
walls were left straight. Consistent with linear theory, the wall boundary con- 
ditions are changed such that the velocity perturbation normal to the wall is 
equal to the slope. The centerline pressure distribution Cp = - 2u/U^ for 
the two cases (one wall or both upper and lower wall modified) are shown in 
Figure (9). It can be noticed that the pressure effect can be felt up to the 
region of the model and can be substantial. Shown in Figure (10) is the lift 
interference 6. When both walls are modified, the lift interference is the 
same (as it should be) as in the case of a straight wall, but when only the 
upper wall is modified the variation of 6 may be significant near the end of 
the slotted region. 


CONCLUDING REMARKS 

The range of validity of the linear representation of the flow fields con- 
sidered here is limited to Mach numbers below the transonic range. Nonlinear 
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FIGURE 8. EFFECT OF INCREASE OF BOTTOM PRESSURE ON OUTFLOW FROM LOWER WALL 
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FIGURE 10. EFFECT OF MODIFICATION OF SOLID WALLS DOWNSTf’EAM OF SLOTS ON INTERFERENCE PARAMETER 




effects are certainly of Importance In transonic Interference, c.f. also 
Reference (10); thus more refined numerical calculations with the non 1 Inear 
transonic equations may be warranted. The present analysis and results, 
although more of a qualitative than a quantitative nature are, however. 

Indicative of important trends of Interference patterns In tunnels with 
slotted test sections of finite length. 

The finite length of the test section appears to be of Importance for 
the definition of interference at the wing location and for the correct in- 
terpretatlon of pitching moment measurements of ccwiplete configurations, when 
the model is placed near the end of the test section, as Is the case with some 
present installations. More significantly, the present results indicate the 
strong sensitivity of the 1 i ‘t Interference to the properties of the slotted 
v#all (a well known fact) and small perturbations of flow conditions which 
could be pre?''Pt in the wind tunnel, such as variations of static pressure 
between lower and upper plenum. These variations can be generated by the loss 
of momentum of high velocity air entering the plenum, by pressure losses of 
air circulating from bottom to upper plenum and, at transonic speeds, by the 
different mass fluxes of air entering the tunnel (above the wing) and entering 
the plenum (below the wing), because of the different temperature of the two 
streams. Also, the admission of stagnant air from the plenum into the tunnel 
decreases the momentum of the tunnel flow and can produce a small pressure rise 
in the flow downstream of the slotted section. This can affect the lift int - 
ference and create axial pressure gradients. 

It is evident that a ^.uant i tat i ve definition of the perturbations induced 
in the tunnel flow by the slotted tert section is needed for a firmer assess- 
ment of interferer.u* effects. These effects can be accounted for and possibly 
corrected, for a given tunnel configuration, if pressure measurements are taken 
during a test at or near the tunnel walls, both along slotted and solid sections. 
These measurements will define directly the wall boundary conditions thus by- 
passing the necessity of using the uncertain homogeneous boundary conditions. 

The present analysis, once further developed to include wing-body configura- 
tions, can permit the determination of interference effects and indicate the 
feasibility of first order corrections by simple modifications of wall charac- 
teristics (for instance, slightly convergent test section, and differential 
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variation of plenum pressure between upper and <ower end of the slotted 
section) when the runnel flow is subsonic. 
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